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Abstract 

' Here we shall find the green's function of the difference equation of 

CO . loop quantum cosmology. To illustrate how to use it, we shall obtain an 

' iterative solution for closed model and evaluate its corresponding Bohmian 

trajectory. 
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O '. 1 Introduction 

^ , In loop quantum cosmology the evolution of the universe has been investigated 

from different points of view. Some of the results are based on the discrete 
quantum domain whose dynamics is described by a difference equation^] E]. 
As one expects the effects of discreteness are most important at small volumes, 
^ , close to the classical singularity. In this limit the evolution of the universe is 

■ completely different from that given by the Wheeler-DeWitt equation. There- 

fore one has to learn how to extract results directly from the difference equation. 
The behaviour of quantum states near the Planck scale is not entirely known 
and only some numerical results exist. In ^ using an isotropic closed model 
with a massless scalar field, the divergent behaviour of the wave function at 
large scale is investigated numerically. One of the most important consequences 
of discrete dynamics is the possibility of evolving from negative values of the 
scale factor eigenvalues to positive ones for all homogeneous models In this 
way the quantum dynamics does not break at the classical singularity. 

At large scales, one can approximate the difference equation with a differ- 
ential equation, where the discreteness leads only to small corrections. This 
is an intermediate semiclassical phase in which the evolution equation takes a 
continuous form There is another analysis called effective classical analy- 
sis. In this approximation one uses the effective matter Hamiltonian in which 
instead of inverse volume operator, it's eigenvalue is used ^IHIEI- In this way 
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one gets the effective classical equation of motion valid at large volumes, de- 
scribing the motion of the position of the wave packet in coordinate time[7]. 
The most recent phenomenological applications of loop quantum cosmology are 
based on these eauationsjH 151 lUl [TU|. Inflationary scenarios without the no- 
graceful exit problem^ |S1 El IHl E|j new perspectives on initial conditions for 
standard inflation|Sl ^2 ^^'^ resolution of the Big crunch problem in closed 
universe^ni are some of them. 

In general the difference equation of loop quantum cosmology has solutions 
which have rapid change in their values and even in their sign when the volume 
changes slightlvfT^. The notion of pre-classicality^^ eliminates such unphys- 
ical solutions which have highly oscillatory behaviour for large scales. The 
physical solutions at large scales should be wave packets which their (WKB 
and/or Bohmian) trajectories are the classical ones. So far this condition on 
the solutions of loop quantum cosmology is done using the generating function 
techniaues fniinilTS] . 

Recentlyjni we have obtained the exact solution of the Hamiltonian differ- 
ence equation for the vacuum case. In this paper, we shall use that vacuum 
solutions, and find the green function of the difference equation. This is an 
important step in obtaining analytical solutions of the difference equation. Be- 
cause using the green function one can obtain at least an iterative solution. We 
shall use this green function to obtain an iterative solution of closed universe in 
loop quantum cosmology in the presence of a massless scalar field. Finally we 
shall investigate the Bohmian trajectories corresponding to this solution, and 
see that these trajectories are compatible with the preclassicality condition. 

2 The Difference equation of cosmology 

Classically, the symmetry of cosmological minisuperspace is introduced by choos- 
ing an isotropic triad = p6f and connection — cS^^. The only physical 
parameters, p and c are related to the more familiar scale factor by the re- 
lations \p\ — and c = ^(^ ~ To) where k is the curvature parameter and 
7 is the Immirizi-Barbero factor. Then, the quantum dynamics of states is 
obtained from the Hamiltonian constraint leading to the following difference 
eauationfTTl lTHllDj : 

in which Vl = 2 — Mqj'^T {T — 1), Cm is the matter {(f>) Hamiltonian, and £p is the 
Planck length. T is the spin connection parameter defined as F = fc/2. ip{(j),£) 
is the coefficient of the expansion of the state in terms of spin network states: 

|V)=^V(0,^)l^> (2) 
e 
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(I£l \ "^^"^ 
j ip are the eigenvalues of the volume operator. Finally \i) 

represents the spin network states for this simple minisuperspace, and it is 
simply the eigenstate of p. 

There is a contribution from the eigenvalues of inverse volume operator in 
the matter Hamiltonian. These are represented by dj^^{€): 

^^i)^dMi) (3) 

where /x G (0, 1) and the positive integer j represent the quantum ambiguities 
arising from constructing quantum operator from the corresponding classical 
expressionism ■ These eigenvalues are bounded at very small scale so the inverse 
volume operator is well defined at classical singularity. For large j the behavior 
of dj_i^{a) shows a peak at a* — {"/ipj/SY^'^ and then approaches to zero 
for all values of /i. Therefore the scale a* determines the size of the scale factor 
below which the geometrical density is significantly different from it's classical 
form. For a << a* the density approaches to zero according to a power law 
behavior dj ^{a) ^ and after passing a peak at a* it decays classically 

If one chooses the matter part to be a massless scalar field (This is the model 
we shall use to demonstrate the use of the green function obtained in the next 
section.): 

C,.. . . 4a-|, (4) 

and assume an oscillatory matter field dependence for the coefficient of the state 
vector to the matter field: 

ip{(t),e)^<S>{£)e"^'^ (5) 

where w is a constant. So the right hand side of the difference equationJQ) can 
be written as: 

_ ^73£3^2(^„(f , </,)V.(0, £) = -eljHyd,,^{e)i:{cly, I) (6) 

where — ^^^^lj^. Using a new variable F{i): 

Fii)^{Ve+,,-Vf,^f,„)e-^'^/'m (7) 
the difference equation can be written as: 

F{£ + Uo) - nF{£) + F{i - Uo) ^ Q{£) (8) 

where: 

Q{£) = -Xfd,,,{£)—^^— (9) 
where A is a dimensionless parameter defined as A = £o^p7'^w^. 
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3 Green function 



In this section we shall find the green function of the difference equation 
By the green function we mean a function satisfying the relation: 

Gij] + A£o) - ngiv) + G{ri - 44) = 5^.» (10) 
The general solution to the equation (jSJ is then can be written as: 

F(^)=^^W+^g(£-f)Q(/) (11) 
I' 

In order to make the above summation meaningful one can assume that G has 
support only on a countable subset of the real line. That is to say, one can 
construct an infinite number of green functions each denoted as Q'^'^°\m) = 
Qim-qo)^ where m is an integer and 770 is a real number. Different values of 770 
leads a green function and a solution covering a subset of the domain of i. 



3.1 Closed universe 

To derive the green function we need the vacuum solutions. From reference [TC] 
we know that for the case of closed universe the vacuum solution is of the form: 

e^^^] with cosh^i - (12) 

4^0 2 ^ ' 

The green function is then can be suggested as (the boundary condition is chosen 
such that the green function goes to zero as \ri\ 00): 

ae^-^ 77 < 

g{^)^\ b 7^ = (13) 

ce-^'' 77 > 

Substituting this in the relation l(Tn|l determines the constants as: 

"^^ = '^= 2sinhUo ^'^^ 

and the that 77, as expected from the discussion above, should comes in steps of 
770 with 770 > Mq. The fact that 770 > Mq would be clear if one substitutes the 
above solution in the relation (|10|) , for 77 < Mq one encounters inconsistency for 
the values of a, b, and c. Therefore 77 could not be smaller than 4£q. It is clear 
that changing 770 in the range [Mq, 8£q\ would cover all the range of |^| > Mq. 
The green function would be simplified then to: 



-I3ria\n—m\ 



n 



(15) 



2sinh4/34 
The solution is thus: 

F(r7^77o) = F(")(77^ryo) + ^ ^'""^ {n~m)Q{f]Qn) (16) 
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As it is noted previously, this solution works only for \l\ > 4^o 
one has to use the difference equation to get: 

F(±e) = Q{Uo ± e) - F{Uq ± e) - VLF{A(.q ± e) 

where < e < 4^o- 

3.2 Flat universe 

For this case the vacuum solution is 

e'"^ and te'°'^ with a = ^ 

where j is an integer number. Again the green function is a linear combination 
of vacuum solutions. Putting the suggestion in the difference equation of green 
function, fixes some of the constants and we are left with the green function: 

where a is some undetermined constant. The appearance of a constant in the 
green function comes from the fact that we have not fix the boundary condition 
in the large i limit. It only demanded that it should be oscillatory. It must 
be noted that again this green function can be used only for |£| > 4£o, and for 
|£| < 4£o one should use the difference equation itself. 



. For 1^1 < 44, 
(17) 



4 Iterative solution 

Using the green functions derived in the previous section, one can obtain an 
iterative solution for the wave function. One can rewrite the equation for F as 
follows in terms of the wave function: 

$(to) = (m) + ^ ^g{n- m)S{n)^{n) (20) 

^ ^ n 

where we have dropped 770 for simplicity and AF(m) — Vm-qo+io ~ ^mrio-to^ 

S{n) = ildj^f,{n) and G = e-''^no{n-m)/2g^ j.^^^ ^^^^ written in terms of 
infinite dimensional matrices as: 

AF$ = Ay$(") + \g-^AV^ (21) 



This has the formal solution: 



AF$ = ^^Ay^^") (22) 
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which can be expanded as: 



= ^^''^ (m) + ^^7^ II - m)5(n)$("' + • • • (23) 

In order to guarantee the convergence of this solution one should choose A as 
small as needed. This makes some limitation on uj. 

In the following we shall apply this general solution to the case of closed and 
flat universe and derive the wave function up to the first order of iteration. 



4.1 Closed universe 

For a closed universe the zeroth order solution is: 

a>L°)(TO) = e--/.^ (24) 

If one ignores the small terms like e~^^^°'^ appearing in the sum needed for 
calculating the first order of iteration, in comparison to 1, and using the different 
behavior of dj,^ in different regimes one can obtain the closed form: 



' A1/(to) \ 2sinh4/34 V^^M 



where 



E(rn) = < , , , V^m* (2 + m)/(2-2m) (26) 

and '^{l^x) — <PhiT{x) / dx^ , in which F is the Gamma function. In the above 
relation m* corresponds to a* . 



4.2 Flat universe 



For this case, the iterative solution is not useful. Because in the summations, 
terms like 1/n appears and this diverges. This clearly does not mean that 
the solution does not exist in this case. The flat isotropic cosmology with a 
massless scalar field investigated by Ashtekar et. al. recently. They used a 
semiclassical state as an initial state and evolve it backward in time and show 
that bigbang is replaced by big bounce. Our iterative method diverges for this 
model, and this means that the above results cannot be obtained from it unless 
one find a way to regularize the iteration. 



5 Trajectories 

Here we shall find the Bohmian trajectories corresponding to the iterative wave 
function of closed universe obtained here. To do this, it is needed to have the 



6 



wave function in the configuration space |c). Using the fact that (c|£) — e'^'^/^, 
one has: 

rio ™ 

Here we choose /(w) = 6{uj — ujo) and g{r]o) = SrjoX (C > 4^o)- Using the 
relation H25|) , one can obtain the wave function in the configuration space up to 
first order. It is plotted in figure 

The corresponding Bohmian trajectory can be found simply by the guidance 
relation. Using relations of reference we have: 

The coordinate time can be eliminated between these equations and a as a 
function of (p derived. The result is plotted in figure 101 ■ Also in this figure the 
classical trajectory is plotted. The classical trajectory can be obtained from the 
classical Hamiltonian constraint P]: 

.2 / 2 „ , , 2 



with — uj one gets the Friedman equation: 

a' + l^^ (31) 

which is independent of 7 as it should be. The solution to this equation is either 

a = constant = (SvrG'tD^/S)"'^/^ or 

t - / , ^ (32) 



which can be expressed in terms of elliptic integrals. 

It is important to note that the scale factor has oscillations of amplitude 
of order of dp. For times far away from the classical singularity, the Bohmian 
trajectory is the classical one up to quantum oscillations (fluctuations) around 
it. This is in accordance with the pre-classicality condition. 

At this end it is fruitful to note some points. First that this is an iterative 
method based on the standard idea of Green's function. Naturally this may be 
of limited use. It is not possible to obtain the exact solution from the iterative 
one, in general. This method is useful! only if the iteration is convergent to the 
exact one. 

Second, which is related to the first, is the reader may ask what is the 
meaning of smallness of A? It is assumed that if one chooses A small enough. 
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which means choosing the matter frequency lo smah enough, it is possible to 
make the iteration convergent. Physically this means that if the matter fre- 
quency is small, the matter contribution to the difference equation would be 
small and one expects that iteration around the vacuum solution is convergent. 
But mathematically, one should examine whether each model is convergent by 
rigour mathematical methods. 

Finally, it must be noted that the method presented here can, in principle, 
be extended to more complicated models. For example, for anisotropic models, 
the only difference is that we have a partial difference equation. The Green's 
function of this equation can be obtained via separation of variables. The result 
would be some kind of multiplication of the isotropic Green's function for each 
direction of the anisotropic model. 

Acknowledgment: This work is supported by a grant from University of 
Tehran. The authors are grateful to Martin Bojowald for useful discusions and 
suggestions. 



References 

M. Bojowald, Phys. Rev. Lett., 86, 5227, 2001, |arXiv:gr-qc/0102069| . 

M. Bojowald, Phys. Rev. Lett., 87, 121301, 2001, la rXiv:gr-qc/0104072| . 

D. Gr een and W. Unruh, Phys. Rev. D, 70, 103502, 2004, 
[arXive |gr^c/0408074 . 

M. Bojowald, Phys. Rev. Lett., 89, 261301, 2002, |arXiv:gr-qc/0206054| . 

M. Bo jowald and K Vandersloot, Phys. Rev. D, 67, 124023, 2003, 
[arXive igr^/03030721 . 

G. Date and G. M . Hossain, Class. Quant. Grav., 21, 4941, 2004, 
|arXiv:gr-qc/040707^ . 

M. Bojowald, P. Singh and A. Skirzewski, Phys. Rev. D, 70, 124022, 2004, 
[arXive:gr-qc/0408094 . 

M. Bojowald, J. E. Lidsey, D. J. Mulryne, P. Singh and R. Tavakol, Phys. 
Rev. D, 70, 043530, 2004, arXiv:gr-qc/0403106 . 

S. Tsujikawa, P. Singh fc R. Maartens, Class. Quant. Grav., 21, 5767, 2004, 
|arXiv:astro-ph/0311015| . 

P. Singh fc A. Toporensky, Phys. Rev. D, 69, 104008, 2004, 

[arXiv:gr-qc/0312110 ; 

G. V. Vereshchagin, JCAP, 0407, 013, 2004, |arXiv:gr-qc/0406108 |. 

M. Bojowald & K. Vandersloot, Invited parallel talk at X*^ Marcel Gross- 
mann meeting, July 20-26, Rio de Jaeiro, 2003, arXiv:gr-qc/0312103 . 



8 



[12] Living Rev. Relativity, 8, 11, 2005. Available online at: 

\protect\vrule widthOpt\protect\href {http : //www . livingreviews . org/lrr-2005-ll}{http : / /wx 

[13] D. Cartin, and G. Khanna, Phys. Rev. Lett., 94, 111302, 2005. 

[14] D. Cartin, G. Khanna, and M. Bojowald, Class. Quant. Grav., 21, 4495, 
2004. 



[15] S. Connors, and G. Khanna, Lanl arXive gr-qc/0509081 



[16] A. Shojai, and F. Shojai, EuroPhys. Lett., 71, 6, 886, 2005. 

[17] A. Ashtekar, and J. Lewandowski, Class. Quant. Grav., 21, R53, 2004. 

[18] A. Ashtekar, M. Bojowald, and J. Lewandowski, Adv. Theor. Math. Phys., 
7, 233, 2003. 

[19] M. Bojowald, H.A. Morales-Tecotl, Lect. Notes Phys., 646, 421, 2004. 
[20] T. Thiemann, Class. Quant. Grav., 15, 1281, 1998. 
[21] M. Bojowald, Pramana J. Phys., 63, 765, 2004, |arXiv:gr-qc/0402053 . 
[22] A. Ashtekar, T. Pawlowski, and P. Singh, Lanl arXive gr-qc/0602086 



9 




i f #( FSrt 

Imaginary Fare 
I •■ • » • i ! < filBrm :Sqiiaf'8J 



Figure 1: The real part, imaginary part and norm squared of the wave function 
of a closed universe in the configuration space and up to first order of it- 
eration. The horizontal axis is Cc/2, and the vertical axis is scaled by the factor: 
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2sinh4/3£o 



^3/2(l-^)_g(^.) ^ ^(1 _ ^(2 + M)/2(l-M) 
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In plotting the above, the oscillatory matter part e*'^'^ is not included. 
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Figure 2: The classical and quantum scale factors as a function of the matter 
field. The horizontal axis is and the vertical axis is Note that 

for times greater than the classicality time (k^^^lo^^"^), the quantum solution is 
oscillations of amplitude ~ 3.Mp with period ~ gT^, where Tp is planck's time. 
For small times, the Bohmian quantum trajectory has no singularity. 



11 



